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TECHNICAL MEMORANDUM 


A CASE STUDY OF ANALYSIS METHODS FOR LARGE DEFLECTIONS 

OF A CANTILEVER BEAM 

INTRODUCTION 


Occasionally the analyst is required to predict the maximum bending deflection and moment of 
cantilever beam where the usual linear methods are insufficient. This occurs in slender beams, such as 
springs, when deflections are large enough that the exact nonlinear Bemoulli-Euler equation must be 
solved. The basic assumption is that the beam is flexible enough to allow large deflections while the 
strains remain elastic. Other assumptions are that the beam is initially straight and of uniform cross 
section and that the bending does not alter the length. This report illustrates analyses via closed-form 
solution and finite-element methods for three common load cases and presents nomograms for easily 
obtained approximate solutions. 


CASE I: Deflection of a Cantilever Beam With Load P, 
Which Remains Vertical, at the Free End 1 


The length of the beam (undeformed or deformed) is given by 

L\ = L—L2 , 


where L2 is the imaginary portion of the beam. 1 Now 


and 


and therefore 


where 


L 2 = J F (P’<Pb) ’ 

L = ±K(p) , 

L, = i {k(p)-f[ P , sin' 1 . 



and E is the material modulus of elasticity and I is the area moment of inertia. Kip) - Legendre s com- 
plete elliptic integral of the first kind and 


Fip4 B ) = J Q 


[dy/lil-p 2 sin 2 y/)2] 


is Legendre’s elliptic integral of the first kind. Where, for this case. 


since 


= sin" 1 (0.7071 Ip) , 


and 


-sin j- 6 B = p { sin <p B ) , 


d B = -±n . 


The value of p must be found by trial and error from the equation for L\. Having found p, the solution 
proceeds as follows: 


h = 


2 p 

k 


the horizontal deflection of the end of the beam is 

S h = L l -h cos <j> B , 


and the vertical deflection is 

8 v = f [E(p4 b ) - \ F{p,<t> B ) - E(p) + lkL] . 
where E(p) is Legendre’s complete elliptic integral of the second kind and 


E{p,<t> B ^=[ (l-p 2 sin 2 y /)2 dy , 


is Legendre’s elliptic integral of the second kind. The geometrical representation of most of the above 
variables is shown in figure 1 . x 
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An example problem (example 2, ref. 1) will now be solved to demonstrate the use of these 
equations. Find the horizontal and vertical deflection of the end of a 150-inch long horizontal cantilever 
of 2 l /2 by 'A inch cross section with a vertical load P = 5 lb applied at the end (E = 30xl0 6 lb/in 2 ). 

/ = btV 12 = (2.5)(0.25) 3 /12 = 3.2552xl0' 3 in 4 , 

El = 9.7656X10 4 lb-in 2 , 

P/El = k 2 = 5.12x10^ in -2 =>k= 1/139.754 in- 1 = 7.1554xl0- 3 in- 1 , 

L, = 139.754 [K(p)-F\p, shr 1 (0.7071/p)]} = 150 in . 

This equation can be solved for p by trial and error in the following manner: 

K{p)-F{p,<p B ) = 139^754 = 1*0733 , 

where 

0 B = sin -1 (0.707 1 Ip) , 

and introducing the expression 

p = sin ^ 6 c (ref. 2) , 

where d c is the rotation at the end of the beam relative to the vertical axis. Initially let d c = 100° which 
is a 10° rotation of the end with respect to the horizontal axis. 

p = 0.766 => <Pb = 67.3766° . 

Use elliptic integral tables to find K(0.766) and F(0. 766, 67. 3766°). 

Kip) = 1.9356 at sin- 1 (0.766) = 50° 
and 

F(p,^b) = 1.3414 at sin- 1 (0.766) = 50° 

(linearly interpolated values between 67° and 68°). 

K{p)-F{p,<pB ) = 0.5942 * 1.0733 , 


therefore try another value of dc- Let 

6 C = 140° =* p = 0.93969 =><pB = 48.8057° . 


Using elliptic integral tables again, 

K(p) = 2.5046 and Fip^s) = 0.95792 , 
K(p)-F(p<t> B ) = 1.54668 * 1.0733 , 
therefore try another value of 6c- Let 
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dc = 120° => p = 0.86603 =><j> B = 54.7348” . 
Then using elliptic integral tables yields 

K(p) = 2.1565 and F(p^ B ) = 1.0783 , 


K(p)-F{p,<t> B ) = 1.0782 ± 1.0733 . 

Now. using these values, h = 242.06 in and L = 139.754(2.1565) = 301.38 in and the horizontal 
deflection is 

S h = 150-242.06 cos (54.7348) = 10.244 in. 

Two additional terms are evaluated from the elliptic integral tables for the vertical deflection 

E(p) = £(0.86603) = 1.21 11 , 

E{p,<p B ) = £(0.86603,54.7348°) = 0.8555 , 

therefore 

& = 2( 1 39.754) { 0.8555-0.5( 1 .0782)— 1 .2 1 1 1+[30 1 .38/(2)( 1 39.754)] }= 51.3 in. 

According to linear theory, the vertical deflection would be 


<S V = 


PL 3 

- = 57.6 in 


and, of course, the horizontal deflection is zero. The explanation for the larger linear deflection is that 
when the load increases gradually from 0 to 5 lb the bending moment is not increasing at the same rate 
since the horizontal movement of the end of the cantilever is toward the fixed end. The linear theory, on 
the other hand, ignores the nonlinear relationship between load and bending moment. This is why the 
principle of superposition is not valid when applying the exact theory. The exact bending moment is 

M max = P(Li-S h ) = 5(150-10.244) = 698.78 in-lb , 


and the maximum bending stress is 


M r 

^max = — = 26,833 lb/in 2 . 

For comparison purposes, this problem was solved using the finite element computer code SPAR. 5 Input 
data and results are contained in appendix A. Maximum deflection, moment, and stress results are 

& = 50.72 in and 4 = 10.72 in . 

Mmax = 696.4 in-lb and tr max = 26,742 lb/in 2 . 

There is less than 5-percent difference in the results of the two methods. 
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CASE II: Deflection of a Cantilever Beam With a Follower Load P, at the Free End, 
Whose Line of Action Remains Perpendicular to the Beam Axis 


The classical analysis of this case rests with the solution of 

El — % - P, cos <p - P 2 sin 0 = 0 , 
ds 2 

where Pi = P sin S B and P 2 = P cos 0 B , <p is the rotation of any point along the beam and 8 B is the angle 
from the horizontal axis of the line of action of P. The classical solution to given below, 

omitting the derivation. 3 The horizontal deflection at any point along the length of the beam is 

§ =1 {cos 0 B [F(p,m)-F(p,n)+2E(p,n)-2E(p,m)]+2p sin 0 B (cos m - cos n)} , 


and the vertical component is 

8 y = j-{2pcosd B (cos m - cos ri) — sin 0 B [F(p,m)-F(p,n)+2E(p,n)-2E(p,m )] } . 
The maximum horizontal deflection at the free end (from the fixed end) is 


x = 1 {cos 6 d [F(p,m)-K(p)+2E(p)-2E(p,m)]+2p sin 6 B cos m} , 
c k 

and the maximum vertical deflection is 

y c = 8 v = i { 2p cos 6 B cos m - sin d B [F(p,m)-K(p)+2E(p)-2E(p,m )] } , 

where 



as in case I. For this case 


where 


p = -sin ^ (8 B +y) , 
8 B +y = - 90° , 


and Y i s the rotation at the free end relative to the horizontal axis. Also 

m - sin -1 [- j sin j 8 B ] , 

and rill 

n = sin 1 1- -p sin (0 B +^)J , 

where $ = y and cos n = 0 at the free end. E(p) = Legendre’s complete elliptic integral of the second 
kind and E(p,m or n) = Legendre’s elliptic integral of the second kind 
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(l-p 2 sin 2 <f>)d</> . 


E{p,m, or n) = 



The solution to this problem can be found by applying the principle of elastic similarity 3 and is 
presented in the following, omitting some derivation. The length of the beam (deformed or undeformed) 
is found to be 

Lx={mp)-F{pjn)} . 

The horizontal and vertical components of the deflection of a point Q(x,y ) along the beam can be found 
in terms of the .yy-coordinate system. The quantities r and q in the x'y '-coordinate system are given as 

r = ^ [E(p,m)-F(p,n)+2E(p,n)-2E(p,m)] , 

where p, m, and n are defined above. At the free end of the beam, point C, 

r c=j [F(p,m)-K(p)+2E(p)-2E(p,m)] , 

q = h{ cos m - cos n) = 2p(cos m - cos n)/k , 


and at the free end 
Hence, at any point Q 


At point C, the free end 
and therefore 


2p 

q c = -f-cosm . 

5 x = r cos 6 B + q sin 8 B , 
Sy = q cos 6 B - r sin 6 B . 


x c = r c cos 0 B + q c sin 0 B , 
S h = E j —x c , 


yc= &M = ( ic cos 0 B~ r C sin 6 B ■ 

The geometrical representation of most of the above variables is shown in figure 2. 

An example problem will now be solved to demonstrate the application of these equations. Find 
the vertical and horizontal deflection at the end of a cantilever beam 20 inches in length with a 0.2 inch 
square cross section and loaded by a force of 5.347 lb at the end which remains perpendicular to the axis 
of the beam. The material modulus of elasticity E of 10x1 0 6 lb/in 2 . 

I = j^ = 1.3333x10^ in 4 , 

9b+Y= -90°, which implies that the load remains perpendicular to the beam axis 
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and 


p = -sin j ( 0 B +y) = 0.7071 , 

k= \fii = 6 - 3327xl0_2in ' 1 * 



Figure 2. Case II variables and geometry. 2 
Now the equation for length L\ may be solved for the angle 6b- 

Li = 15.791 [K{p)-F{p,m)] = 20 . 

K{p) = K (sin- 1 p) = £(45°) = 1.8541 from complete elliptic integral table 
F(p,m) = 1.8541 - = 0.587556 . 

This relationship is satisfied if m = 32.765° from tables of elliptic integrals of the first kind. Therefore, 
knowing p and m, from the elliptic integral tables 

E(p,m) = 0.55692 , 

Eip) = 1.3506 , 

and from the equation for m, d B = -45° which means that the rotation of y of the free end is equal to 
-45°. Now all the terms are available to evaluate r c and q c 

r c = 5.066, 
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q c = 18.77885 . 

Using the absolute value of 6g, x a and y c are calculated to be 


jc c 5.066(0.7071)+ 18.77885(0.7071) = 16.861 inches 
and 

8 h = 20.0-16.861 = 3.139 inches 
y c =8 v = 18.77885(0.707 1)— 5.066(0.7071) = 9.696 inches . 

The bending moment A/ max at the fixed end is equal to Pq c . 

M m « = 5.347(18.77885) = 100.41 in-lb . 

This problem was solved using the finite element computer code SPAR. 5 A follower force was applied 
at the free end of the beam. Input data and results are contained in appendix B. Maximum deflection, 
bending moment, and stress results are 

& = 9.697 inches and 5^ = 3.137 inches , 

M max = 100.42 in-lb and (W = 75,310.0 lb/itf. 

There is no significant difference between the results of the numerical and exact theoretical methods for 
this example. 


CASE III: Deflection of a Cantilever Beam With Uniformly Distributed Load 

Whose Direction Remains Vertical 4 


In this case, as in the previous cases, only beams of uniform cross section are considered. The 
basic assumptions are that the deformations are elastic and that the bending does not alter the length of 
the beam. 

Cases I and II obtain results for a concentrated load. In the case of a uniformly distributed load, 
the Bernoulli-Euler equation, 

J_ _ dd_ _ M_ 

P ds El * 

gives, with the origin at the free end of the beam, 

4¥- = -ws cos 6= El , 
ds ds 2 

where p is the deflected beam curvature, 8 is the rotation, s is the distance to any point along the length, 
M is the bending moment, E is the modulus of elasticity, I is the cross sectional moment of inertia, and w 
is the load per unit length. 

This equation does not lend itself to any simple solution, therefore let 
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0=Zfl/;y=i na ^ ; y-y = £ n(n-l)<v" 2 . 

n = 0 05 n = 1 ds " = 2 

The boundary conditions s = 0,8= a, dO/ds = 0 give oq = or ; 01 = 0. Then 


(cos a cosT - sin a sin T) , T = £ a^s" , 
ds 2 EI "= 2 


where or is the rotation of the beam at the free end. 


Expand cos T and sin T as 


- T 2 T 4 T_ 

cos T - 1 - 2j + 4| 5| + ••• ’ 

t'3 r p ^ t7 

sin T = T- -jp + - 7-7 ~~j\ + * 

and substitute into the previous equation. Equating the two expressions for — -y and substituting the 

above expression for T, the coefficients a n may be determined by equating like coefficients of s. The 
coefficient o 2 of the constant term is equal to zero. Coefficients of 5 through s 1 are: 

a 3 = - U ~ > (coefficient of 5) 

a 4 = a 5 = 0 , (coefficients of s 2 and s 3 ) 


_ a 3 w sin or 
° 6 = 30£7 


— s * n — c ° s — , (coefficient of 5 4 ) 
180E 2 / 2 


a n = a % = 0 , (coefficients of 5 5 and s 6 ) 


and 

Thus, 

Since dy/ds = sin 6, 


_ a 6 w sin or 
~ TIE! 


(coefficient of s 1 ) 


a 3n+4 = a 3n+5 = 0 f° r n = 0, 1, 2, 3, ... . (ref. 4) 
6 = a + £ a 3k s 3k . 

* = i 

y = f sin 0 J5 = Tj sin or + T, cos or , 

Jo 

r _ 

14 10 

^ a 3 s 4 ' a tf 1 ( 09 - 03 / 6 ) 5 10 

2 2~ — + 7 + 10 
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Similarly, 


jc = r, cos a - T 2 sin a . 

Finally, solving for M in the Bemoulli-Euler equation and differentiating the expression for 8 
with respect to s gives 

M = El £ 3ka 3k s 3k ~ l , 
a 3k L 3k , 

where, for this case, the imaginary beam length is zero and L is the actual beam length. With s = L,a = 
wL 3 /6El, the first approximation to the maximum deflection is 

1 r 4 

y ± s sin a + ^ a 3 s 4 cos a - , (sin a - a and cos a ± 1) 

and the moment is 

M max ±EI(3a 3 s 2 )± , 

which are in agreement with small deflection theory. 

The following example illustrates these equations. Given a thin cantilever beam 20 inches long 
and constant cross section of 0.2 inches square, find the uniformly distributed load which causes a 1.0 
radian rotation at the free end and find the maximum deflection and stress for that load. The modulus of 
elasticity is 10xl0 6 lb/in 2 . 

/ = lT = 1 - 333xl0 ^ in4 • 

is the area moment of inertia and 

a = - £ a 3k L 3k = -a 3 L 3 -a 6 L 6 -a 9 L 9 ... 

_ wT 3 cos a + w 2 L 6 cos a sin a w 3 L 9 cos a sin 2 a 
180£ 2 / 2 12,960 £ 3 / 3 

Let C = —£f an d substitute in the above equation. After moving a to the right hand side and dividing by 
"12:960 ~ ex P resS30n 

=0. 

sm u sin a sin a cos a 

Solving this equation for £ with a = 1.0 radian gives £ = 8.7422. Therefore, the uniformly distributed 
load w is equal to 1.45667 lb/in. Continuing with the solution 
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a 3 =-9.84xl0- 5 , 

a 6 =-3.01623x10-9 , 

a 9 =— 3.852139xl0- 14 , 
and 

T x = 20.0-0.88526-0.30392 = 18.81082 
T 2 =-3.936-0.55154+0.12316 = ^1.36438 . 

Hence, the maximum vertical deflection is 

y L = 8y= 18.81082 sin (1.0) + (-4.36438) cos (1.0) = 13.47 in. 

The horizontal distance of the free end from the fixed end is 

x L = 18.81082 cos (1.0) -(-4.36438) sin (1.0)= 13.836 in. 

Therefore, the maximum horizontal deflection is 

8 h = 20.0-13.836 = 6.164 in. 

The bending moment at the fixed end is 

M max = EI(3a 3 L 2 +6a 6 L 5 +9a 9 L*+...) 

= 1.333x10 s (-0.11808-0.05791-0.00888) 


= -246.43 in-lb . 

Thus, the maximum bending stress is 

^max = = 184,822.5 lb/in 2 . 

The SPAR 5 finite element results are: 

8 V = 13.314 in , 

8 h = 6.0485 in , 
and 

M max = -227.27 in-lb . 

The difference in sign is due to the different coordinate system origins. Note that the finite element 
results are all less than the closed form solution with the largest difference being in A/ max (± 8 percent). 
This difference may be reduced if the total load (wL = 29.13 lb) is concentrated at midlength and 
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multiplied by the deformed distance to the fixed end (8.225 in) to give M max = -239.62 in-lb. Complete 
results are contained in appendix C. 

These may be compared with linear beam theory results: 

<5, = 21.8555 in, 

8h = 0.0 in , 
and 

M max = 291.33 in-lb. 


SUMMARY 


Solutions obtained by the methods presented in this paper are accurate, but lengthy and time con- 
suming. Therefore, the following nomograms, for cases I, II, and HI, are presented in figures 3, 4, and 5, 
respectively, for the analyst who requires a faster but less accurate result. They are derived from the 
closed-form solution equations. Abscissa (x-axis) and ordinate (y-axis) values are nondimensional. 

Abscissa parameters are noted in the legend to the right of each nomogram. Ordinate values are for 

o El 

wL* 

figures 3 and 4 and for figure 5. The parameter L for both axes is equal to the actual beam length 
L\. See the list of symbols for parameter definitions. 

Finite element method results are nearly identical to the closed-form solutions as shown in all 
three cases. This instills confidence in each method. 

Other load cases for which closed-form solutions were found are thermal, 6 two concentrated 
vertical loads, 3 and end moment plus vertical force. 1 
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Figure 3. Cantilever beam with vertical end force. 
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Figure 5. Cantilever beam with uniformly distributed load. 
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APPENDIX A 


Case I Finite Element Model Input and Results 


*XQT TAB 
START 51 2 4 6$ 

TITLE 'LARGE DEFLECTION OF CANTILEVER BEAM WITH VERTICAL END FORCE 
JLOC$ 

1 0. 0. 0. 0. 0. 150. 51 1$ 

MATC$ 

1 30. +6 .3 .298 1.0-5$ 

MREF$ 

11111.$ 

BA$ 

DSY 1 3.2552-1 0. 3.2552-3 0. .625 1 . 302-2$ . 25X2 . 5 RECT 
1. 1. 1. .125 1.25 .125 -1.25 -.125 1.25 -.125 -1.25$ 

CONS=l$ 

ZERO 1 3 5:1$ 

*XQT ELD 
E21$ 

121 50$ 

*XQT E 
*XQT EKS 
*XQT U1 
*(GNA OPTIONS) 

!SS=. 25$ INITIAL STEP SIZE 
! SSMIN= . 1$MINUMUM STEP SIZE 
*XQT AUS 

ALPHA: 1 CASE TITL 3$ 

l'VERTICLE FORCE REMAINS VERTICLE AT JOINT 51,WL»*2/EI=1 . 152 
$1' FORCE FOLLOWS ROTATION OF JOINT 51, (WL**2/EI) =1 . 604 
SYSVEC : 2 APPL FORC : 1=1 : J=51 : -5 . 000$ <WL**2/EI) =1 . 152 
$SYSVEC : 2 POLL FORC : 1=1 : J=51 : -5 . 3470$ (WL**2/EI) =1. 604 
$1 ' 5LB FORCE APPLIED AT 45DEG TO BEAM 

$SYSVEC:2 APPL FORC : 1=1 : J=51 : -3 . 5355 : 1=3 : J=51 : -3 . 5355$ 

•CALL (18 NL JCL) 

•CALL (GNA MAIN) $ PERFORM NONLINEAR ANALYSIS 
•CALL (GNA MOTI ) $ COMPUTE NONLINEAR DISPLACEMENTS 
•CALL (GNA REAC ) $ COMPUTE NONLINEAR REACTIONS 
•CALL (GNA PCH) $ PRINT CONVERGENCE HISTORY 
*XQT VPRT 
FORMAT =4$ 

TPRINT NONL MOTI 31$ 

TPRINT NONL REAC 3 1$ 

*XQT DCU 
TOC 1$ 

TOC 2$ 

*XQT EXIT 
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1NONL MOTI 

VERTICAL FORCE REMAINS VERTICAL AT JOINT 51,WL 


0 JOINT, 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


1 

O.OOOOOOOOOOOE+OO* 
-0. 3185775468 3E -01 
-0 . 12649 6 8723 0E+ 00 
-0 . 2 8 24 9 6 15 2 3 6E+ 00 
-0 . 49841337723E+00 
-0.77278915413E+00 
-0 . 110415 052 18E+01 
-0 . 14910143236E+01 
-0.1931890 34 96E+01 
-0.24 2528425 14 E+01 
-0.29697002 34 6E+01 
-0.35636435348E+01 
-0.42056 22 68 14E+01 
-0.48941515589E+01 
-0.562775 12 69 6E+01 
-0.64049518078E+01 
-0.72 24 293551 1E+01 
-0.80843285782E+01 
-0 . 89 83 621820 6E+01 
-0 . 99207520554E+01 
-0.10894312749E+02 
-0.11902912758E+02 
-0 . 12945176894E+02 
-0 . 14 0197463 6 8E+02 
-0 . 15 12 527910 5E-f 02 
-0.16260 44 995 6E+02 
-0 . 17423950802E+02 
-0.186 14 49055 7E+02 
-0.1983 07 9508 3E+02 
-0 . 21071607019E+02 
-0. 2233568552 9E+02 
-0. 2362180597 9E +02 
-0 . 24928759 544E+0 2 
-0.2625535 27 52E+02 
-0 . 27600406971E+02 
-0 . 28962757845E+02 
-0 . 30341254675E+02 
-0 . 31734759763E+02 
-0 . 33142147708E+02 
-0 . 34562304672E+02 
-0 . 35994127607E+02 
-0 . 37436523454E+02 
-0 . 38888408318E+02 
-0.40348706615E+02 
-0.41816350 19 8E+ 02 
-0.43290277467E+02 
-0.44769432463E+02 
-0.46252763949E+02 
-0 . 4773922447 5E+0 2 
-0 . 49227769441E+02 
-0 . 50717356151E+02 


3 

O.OOOOOOOOOOOE+OO* 
-0.16914902874E-03 
-0. 1662255783 3E-02 
-0 . 572092220 8 9E-02 
-0 . 13500994217E-01 
-0.26074280 12 9E-01 
-0.44430409009E-01 
-0.6 94 7880976 9E-01 
-0 . 10205079293E+00 
-0.14 29017180 9E+00 
-0 . 19271323119E+00 
-0 . 25209555687E+00 
-0 . 32158983217E+00 
-0 .4016 704691 5E+00 
-0.49274753477E+00 
-0 . 59516913773E+00 
-0 . 70922381277E+00 
-0 . 83514289392E+00 
-0 . 97310286921E+00 
-0 . 11232277 10 1E+01 
-0.12855911700E+01 
-0.14 60219 04 65 E+01 
-0 . 164709 14 032E+01 
-0 . 184 614474 8 0E+01 
-0 . 20 5727416 3 6E+01 
-0.22803353896E+01 
-0 . 251514685 30E+01 
-0 . 276149164 64E+01 
-0 . 30191194532E+01 
-0 . 32877484193E+01 
-0 . 35670669711E+01 
-0 . 38567355808E+01 
-0.41563884789E+01 
-0 . 44656353159E+01 
-0 . 47 84062773 6E+01 
-0 . 5 11 12 3 6 12 8 2E+ 01 
-0 • 54467007664E+01 
-0 . 57899836569E+01 
-0 . 61405947796E+01 
-0 . 649 80285 13 6E+01 
-0 . 68617 649 883E+01 
-0 • 72312713982E+01 
-0.76060032856E+01 
-0 . 79854057924E+01 
-0 . 83689148861E+01 
-0 . 87559585596E+01 
-0 . 914595 80 111E+01 
-0 . 95383288050E+01 
-0 . 99 32482 0168E+01 
-0 . 10327825367E+02 
-0 • 107237 64343E+02 


ID= 3/ 

* *2/EI=l . 152 
5 

0.00000000000E+00* 
-0. 2116 24 9696 7E -01 
-0.41864320417E-01 
-0.62 10576906 6E -01 
-0 . 818 873257 54E- 01 
-0 . 10 12 09 6444 1E+ 00 
-0 . 12007353726E+00 
-0 . 13 847 9 9 624 2E+ 00 
-0 . 15 6430011 65E+ 00 
-0 . 17392489 8 60E+ 00 
-0 . 19 09 6 5 94 7 2 7E+ 00 
-0 . 20755458086E+00 
-0 . 22369231102E+00 
-0 . 23938072739E+00 
-0 • 25462148752E+00 
-0 . 26941630721E+00 
-0 . 283766 95 115E+00 
-0 . 2976752239 3E+00 
-0 • 31114296150E+00 
-0 . 3241720228 6E+00 
-0 • 3 3 67 64 2 8 22 0E +00 
-0 • 34892162137E+00 
-0 . 360645 9226 5E+00 
-0 . 3719 39 06 19 1E+ 00 
-0 . 38280290212E+00 
-0 . 39323928707E+00 
-0. 4032500355 5E +00 
-0 . 41283693576E+00 
-0 . 42200173997E+00 
-0.43074615963E+00 
-0 . 43907186059E+00 
-0.44698045870E+00 
-0.45447351563E+00 
-0 .46155253 50 0E+00 
-0 .4682189 587 0E+00 
-0 .474474163 5 0E+ 00 
-0.48031945784E+00 
-0.48575607894E+00 
-0.49078519006E+00 
-0.49540787799E+00 
-0.49962515075E+00 
-0 . 50343793553E+00 
-0 . 50 6 847 07 6 7 8E+ 00 
-0 . 5098533 34 52E+00 
-0 . 5 1245 7 3 82 8 3E+ 00 
-0 • 5146598085 8E+00 
-0 . 5 164 6 1110 2 0E+ 00 
-0 . 51786169 684E+00 
-0 . 51886188750E+00 
-0 . 51946 19 104 6E+00 
-0 . 51966190289E+00 


1 / 1 
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1NONL REAC ID _ 3 / 

VERTICLE FORCE REMAINS VERTICLE AT JOINT 51, WL**2/EI=1 . 152 
0 JOINT, 135 

1 0. 5000000000 3E+01* -0 . 25268581005E-07* 0 . 69638117838E+03* 

2 0.366908 92 84 IE-08 -0 . 74119820420E-07 0 . 14224497136E-08 

3 0.64273665857E-08 -0 . 59487510396E-07 0 .45147316996E-08 

4 -0.83442991160E-08 0 . 18679165254E-06 0 . 93059497885E-08 

5 -0.110 67 002853E-07 0 . 14089676045E-06 0 . 15839759726E-07 

6 0. 4396098302 3E-08 -0 .47157846861E-07 0 .22981112124E-07 

7 0.460797 67 976E-07 -0 . 34994293571E-06 0 . 30788214644E-07 

8 -0 . 46776936667E-07 0 . 37285039037E-06 0 . 41447492549E-07 

9 0 .40927636487E-08 -0 . 19209443135E-07 0 . 52939867601E-07 

10 0. 6199823465 9E-07 -0 . 33212069083E-06 0 . 62398612499E-07 

11 0.28959021320E-06 -0 . 13972768631E-05 0 . 71911927080E-07 

12 -0.45219339283E-06 0 . 22184180677E-05 0 . 86478394223E-07 

13 0 . 22711591727E-06 -0 . 99310839715E-06 0 . 10020448826E-06 

14 0.9704839560 IE-07 -0 . 34056570586E-06 0 . 11163683666E-06 

15 -0 .74678942571E-06 0 . 28669526852E-05 0 . 13378121366E-06 

16 0 . 56335676347E-06 -0 .20977331591E-05 0 . 15072146198E-06 

17 -0.40877 17 6819E-06 0 . 13891253109E-05 0 . 16328704078E-06 

18 0. 6642468618 6E-06 -0 . 21892571347E-05 0 . 17038837541E-06 

19 -0.23180 173345E-06 0 . 76486211739E-06 0 . 17915408534E-06 

20 -0.151865008 04E-06 0 . 45445375184E-06 0 . 19668186724E-06 

21 0.27901017496E-06 -0 . 78814326746E-06 0 . 20764309738E-06 

22 0. 766350 97407E-06 -0 . 20648057751E-05 0 .20850166038E-06 

23 -0 . 17863103046E-05 0 . 47841309528E-05 0 . 22809763323E-06 

24 0 . 31924620270E-06 -0 . 87805593201E-06 0 . 25432018447E-06 

25 0 . 16538629210E-05 -0 .41145539368E-05 0 .24293876777E-06 

26 -0 . 12766291317E-05 0 . 31440978928E-05 0 . 24980181479E-06 

27 -0. 3579335500 6E-06 0 . 81169097785E-06 0 . 26901398087E-06 

28 0.14 6150 12 557E-05 -0 . 33550627957E-05 0 .25135341275E-06 

29 -0 . 18538326464E-05 0 . 41735536966E-05 0 .26310954127E-06 

30 0.2357258279 IE-05 -0 . 51586649736E-05 0 . 26704765332E-06 

31 -0 . 28504115499E-05 0 . 61073636789E-05 0 .26730049285E-06 

32 0.21206145893E-05 -0 .44715303057E-05 0 . 27229089028E-06 

33 -0.64698108018E-06 0 . 13434514541E-05 0 . 25650751923E-06 

34 0.55747 92 873 IE-06 -0 . 11334025390E-05 0 .24536075216E-06 

35 -0. 109692740 17E- 05 0 . 21946723913E-05 0 .24617384042E-06 

36 -0.32130 92 699 6E-06 0 . 61182657252E-06 0 .26677662390E-06 

37 0.68545 54 366 3E-06 -0 . 13588690658E-05 0 .24983910407E-06 

38 0.150389 94 88 5E -05 -0 . 28481178322E-05 0 . 21968662622E-06 

39 -0 . 15744060593E-05 0 . 29410096701E-05 0 . 19704930310E-06 

40 0 . 60735658795E-06 -0 . 10983124328E-05 0 . 20132392820E-06 

41 -0 . 11038678933E-05 0 . 19922122131E-05 0 . 19832623366E-06 

42 0. 173590 1472 6E-05 -0 . 31660113165E-05 0 . 16631474864E-06 

43 -0.1988307 14 5 IE-05 0 . 35696666179E-05 0 . 14507349988E-06 

44 0 .20173751072E-05 -0 . 36301386149E-05 0 . 11937481759E-06 

45 -0 .20609390727E-05 0 . 36645720160E-05 0 . 10267331163E-06 

46 0 .22568194052E-05 -0 . 40014328361E-05 0 . 94547885965E-07 

47 -0.1660 116973 IE-05 0 . 29174758918E-05 0 .79199708125E-07 

48 -0.15050 84 5840E-05 0 . 26605719123E-05 0 . 90053390522E-07 

49 0 .43321427263E-05 -0 . 77020095101E-05 0 . 21684513740E-07 

50 -0.32892274111E-05 0 . 58249666746E-05 -0 . 11591055227E-07 

51 -0 .49999991927E+01 -0 . 14412546392E-05 0 . 18818013814E-07 
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APPENDIX B 


Case II Finite Element Model Input and Results 


*XQT TAB 
START 21 2 4 6$ 

TITLE 'LARGE DEFLECTION OF CANTILEVER BEAM WITH FOLLOWER END FORCE 
JLOC$ 

1 0 . 0 . 0 . 0 . 0 . 20 . 21 1 $ 

MATC$ 

1 10. +6 .3 .101 1.0-5$ 

MREF$ 

11111.$ 

BA$ 

GIVN 1 1.333-4 0. 1.333-4 0. .04 2 . 256-4$ . 2X.2 RECT 
CONS=l$ 

ZERO 1 3 5:1$ 

•XQT ELD 
E21$ 

t 1 2 1 20 $ 

•XQT E 
•XQT EKS 
•XQT ADS 

ALPHA :1 CASE TITL 4$ 

1 ' FORCE FOLLOWS ROTATION OF JOINT 21, (WL**2/EI) =1 . 604 
SYSVEC : 2 FOLL FORC : 1=1 : J=21 : -5 . 3470$ (WL*»2/EI) =1 . 604 
•XQT VPRT 
FORMAT=4$ 

TPRINT 2 FOLL FORC 1$ 

•XQT U1 

* (GNA OPTIONS )$ 

!SS=. 25$ INITIAL STEP SIZE 
I SSMINs . 1$MINDMUM STEP SIZE 
•CALL (18 NL JCL)$ 

•CALL (GNA MAIN) $ PERFORM NONLINEAR ANALYSIS 
•CALL (GNA MOTI ) $COMPUTE NONLINEAR DISPLACEMENTS 
•CALL (GNA REAC ) $COMPUTE NONLINEAR REACTIONS 
•CALL (GNA PCH) $ PRINT CONVERGENCE HISTORY 
•XQT VPRT 
PORMAT=4$ 

TPRINT NONL MOTI 4 1$ 

TPRINT NONL REAC 4 1$ 

•XQT DCU 
TOC 1$ 

TOC 2$ 

•XQT EXIT 
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1NONL MOTI ID* 4/1/ 1 

FORCE FOLLOWS ROTATION OF JOINT 21, (WL**2/EI) =1 . 604 
0 JOINT, 135 

1 0 . 00000000000E+00* O.OOOOOOOOOOOE+OO* O.OOOOOOOOOOOE+OO* 

2 -0 . 37157 97349 8E- 01 -0 . 69969773255E-03 -0 .73860944922E-01 

3 -0 . 14646568 10 9E+00 -0 . 67001513649E-02 -0 . 14468755475E+00 

4 -0 . 32427 657 011E+00 -0 . 22643220231E-01 -0 .21230195271E+00 

5 -0 . 56654657690E+00 -0 . 52441447316E-01 -0.27654930043E+00 

6 -0 . 86894995208E+00 -0 . 99267856809E-01 -0 . 33729667836E+00 

7 -0.12269858297E+01 -0 . 16556585669E+00 -0.39443165037E+00 

8 -0.16360727156E+01 -0 . 25307563728E+00 -0 .44786059467E+00 

9 -0 . 20916291373E+01 -0 . 36287334164E+00 -0 .49750687929E+00 

10 -0 . 25891396451E+01 -0 .49541942456E+00 -0 . 54330895366E+00 

11 -0 . 31242061404E+01 -0 . 65061294532E+00 -0 . 58521841939E+00 

12 -0 . 36925851152E+01 -0 . 82784898817E+00 -0 . 62319813381E+00 

13 -0 . 42902118230E+01 -0 . 10260769158E+01 -0 . 65722039046E+00 

14 -0 . 49132126763E+01 -0 . 12438576857E+01 -0 . 68726521161E+00 

15 -0 . 55579073104E+01 -0 . 14794189620E+01 -0.71331877963E+00 

16 -0 . 62208017952E+01 -0 . 17307072172E+01 -0.73537202717E+00 

17 -0 . 68985744385E+01 -0 . 19954364139E+01 -0 .75341940002E+00 

18 -0 . 75 88055 5437E+01 -0 . 22711332009E+01 -0 .76745780185E+00 

19 -0.82862023713E+01 -0 . 25551788255E+01 -0 .77748572680E+00 

20 -0.89900704394E+01 -0 .28448481858E+01 -0 .78350258303E+00 

21 -0 . 96967821874E+01 -0 . 31373466046E+01 -0 .78550820899E+00 

1NONL REAC ID* 4/1/ 1 

FORCE FOLLOWS ROTATION OF JOINT 21, (WL**2/EI) *1.604 
0 JOINT, 135 

1 0 . 37805317898E+01* 0 . 37812848421E+01* 0 . 10041585312E+03* 

2 -0.3605 133903 0E-05 0 . 80082172810E-05 0 . 10538043625E-04 

3 -0 . 36758077268E-05 0 . 57586682639E-05 0 . 84672487901E-05 

4 -0 . 35562574175E-05 0 .40763069665E-05 0 . 69430980147E-05 

5 -0 . 33287210750E-05 0 .27880235418E-05 0 .58234495555E-05 

6 -0.30727364155E-05 0 . 18756233026E-05 0 .49970376494E-05 

7 -0 . 277 9 12 58 14 IE- 05 0 . 11335455356E-05 0 .43781419663E-05 

8 -0 . 25004520902E-05 0 . 60761937755E-06 0 . 39015853872E-05 

9 -0 . 22134 82 0304E-05 0 . 18376311476E-06 0 .35194223074E-05 

10 -0 . 19275 1875 15E- 05 -0 .15286477211E-06 0 . 31957883948E-05 

11 -0. 1678023524 IE-05 -0 . 36599938256E-06 0.29054508559E-05 

12 -0 . 14201456468E-05 -0 .55291329558E-06 0 .26305763186E-05 

13 -0 . 11698341180E-05 -0 . 69208105134E-06 0 .23588468139E-05 

14 -0 . 99904874661E-06 -0 .70954232548E-06 0 .20832990231E-05 

15 -0 . 77717252442E-06 -0 .78031921873E-06 0.17983695670E-05 

16 -0. 645 16 82989 5E-06 -0 .73406904074E-06 0 . 15023351807E-05 

17 -0 . 50128378716E-06 -0 .70083995827E-06 0 . 11957694142E-05 

18 -0.44193721944E-06 -0 . 57964617921E-06 0 . 87838986929E-06 

19 -0. 3 129 040385 9E-06 -0 . 52649294086E-06 0 . 55287887335E-06 

20 -0. 3115456809 3E-06 -0 . 34264307622E-06 0.22326662474E-06 

21 -0 . 37804968735E+01 -0 .37813031365E+01 0 .11187353266E-08 
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APPENDIX C 

Case III Finite Element Model Input and Results 


platan*** p*gk blank not film ft 


PAGE 


.INTENTIONALLY BLANK 
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•XQT TAB 
START 21 2 4 6$ 

TITLE 'VERIFY CANT. BEAM WITH UNIFORM DIST. LOAD LARGE DISP. 

JLOC$ 

1 0 . 0 . 0 . 0 . 0 . 20 . 21 1 $ 

MATC$ 

1 10.+6 .3 .101 1.0-5$ 

MREF$ 

11111.$ 

BA$ 

GIVN 1 1.333-4 0. 1.333-4 0. .04 2 . 256-4$ . 2X. 2 RECT 
CONS-l$ 

ZERO 1 3 5:1$ 

•XQT ELD 
E21$ 

1 2 1 20$ 

*XQT E 
*XQT EKS 
•XQT AUS 

ALPHA: 1 CASE TITL 4$ 

1 ' UNIFORMLY DISTRIBUTED LOAD, (WL* *3 /El) -8 .7422$ 

SYSVEC : 2 APPL FORC : 1-1 :J»1 : 21 :-. 728335 :J=2, 20 : -1.45667$ (WL**3/EI) -8 . 7422 
$SYSVEC : 2 APPL FQRC: 1-1 : J-ll : -29 . 1334$ (WL**3/EI) =8 .7422 
*XQT VPRT 
FORMAT-4 $ 

TPRINT 2 APPL FORC 1$ 

*XQT U1 

* (GNA OPTIONS )$ 

I SS- .25$INITIAL STEP SIZE 
I SSMIN- . 1$MINUMUM STEP SIZE 
*CALL (18 NL JCL) $ 

•CALL (GNA MAIN) $ PERFORM NONLINEAR ANALYSIS 
•CALL (GNA MOTI ) $COMPUTE NONLINEAR DISPLACEMENTS 
•CALL (GNA RBAC ) $COMPUTE NONLINEAR REACTIONS 
•CALL (GNA PCH) $ PRINT CONVERGENCE HISTORY 
•XQT VPRT 
FORMAT-4 $ 

TPRINT NONL MOTI 4 1$ 

TPRINT NONL REAC 4 1$ 

•XQT DCU 
TOC 1$ 

TOC 2$ 

•XQT EXIT 
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1NONL MOTI ID* 4/ 1/ 

UNIFORMLY DISTRIBUTED LOAD, (WL**3/EI) »8 .7422$ 

0 JOINT, 135 

1 0 . OOOOOOOOOOOE+OO* 0 .00000000 OOOE+OO* 0 . OOOOOOOOOOOE+OO* 

2 -0.81538520445E-01 -0 . 33239993945E-02 -0.15987638662E+00 

3 -0 . 31066211238E+00 -0 .29910802471E-01 -0 .29929428526E+00 

4 -0 . 66394461155E+00 -0 . 94370047518E-01 -0 .41992736056E+00 

5 -0.11195835966E+01 -0 .20417463579E+00 -0 .52358954518E+00 

€ -0 . 16584046886E+01 -0 . 36171831143E+00 -0 . 61209148738E+00 

7 -0 . 22641411169E+01 -0 . 56601274444E+00 -0.68715456592E+00 

8 -0 . 29233073868E+01 -0 . 81397242615E+00 -0.75036631988E+00 

9 -0.36248825814E+01 -0 . 11013321818E+01 -0.80316359804E+00 

10 -0.43599328755E+01 -0 . 14232741968E+01 -0 . 84683307687E+00 

11 -0.51212429944E+01 -0 . 17748413757E+01 -0 . 88252187638E+00 

12 -0.59029891089E+01 -0 . 21512011733E+01 -0.91125342951E+00 

13 -0 . 67004644591E+01 -0 .25478087340E+01 -0 . 93394551116E+00 

14 -0 . 75098579908E+01 -0 . 29605044488E+01 -0 . 95142852387E+00 

15 -0.83280812514E+01 -0 . 33855701642E+01 -0 . 96446291227E+00 

16 -0.91526369711E+01 -0 . 38197602179E+01 -0 .97375506377E+00 

17 -0.99815225441E+01 -0 .42603177514E+01 -0 .97997134508E+00 

18 -0.10813162069E+02 -0 .47049828067E+01 -0.98375009348E+00 

19 -0 . 11646361236E+02 -0 . 51519960813E+01 -0 .98571147435E+00 

20 -0 . 12480279932E+02 -0 . 56001005123E+01 -0 .98646516492E+00 

21 -0.13314417842E+02 -0 . 6048541824 7E+01 -0.98661584764E+00 

1NONL REAC ID- 4/ 1/ 

UNIFORMLY DISTRIBUTED LOAD, (WL**3/EI) »8 .7422$ 

0 JOINT, 135 

1 0 . 28405064739E+02* 0 .37493751393E-06* 0 .22727046458E+03* 

2 -0.14566690868E+01 0 .40573603144E-07 0 . 24162636691E-06 

3 -0 . 14566693891E+01 0 . 17144086151E-06 0 . 10343565009E-05 

4 -0 . 14566696583E+01 0 . 18945757674E-06 3 . 15467057892E-05 

5 -0 . 14566698756E+01 0 . 17425380250E-06 0 . 18266136976E-05 

6 -0 . 145667000 5 5E+01 0 . 92202856905E-07 0.19306394279E-05 

7 -0 . 14566701042E+01 0.40793789680E-07 0.19090766727B-05 

8 -0 . 14566701539E+01 -0.23910054381E-07 0.18009495761E-05 

9 -0 . 14566701822E+01 -0.70063922664E-07 0 . 16351889371E-05 

10 -0 . 14 56670184 9E+01 -0 . 10873932155E-06 0 . 14357237887E-05 

11 -0 . 14566701863E+01 -0.12625688228E-06 0 . 12203508959E-05 

12 -0 . 14566701684E+01 -0 . 14065609122E-06 0 . 10027918051E-05 

13 -0.14566701445E+01 -0 . 14568882634E-06 0 .79395965713B-06 

14 -0 . 14566701346E+01 -0 . 12832246199E-06 0 . 60264358126E-06 

15 -0.14566701210E+01 -0 . 11230199704E-06 0 .43350303258E-06 

16 -0 . 14566701227E+01 -0 .75494916985E-07 0.29091620490E-06 

17 -0 .14566700643E+01 -0 . 82225238174E-07 0 . 17748456571E-06 

18 -0 . 14566700444E+01 -0 . 60133000801E-07 0.92594305556E-07 

19 -0 . 14566700701E+01 -0 . 81951867514E-08 0.36932846115E-07 

20 -0 . 14566700385E+01 0.66707173119E-09 0 . 81222477633B-08 

21 -0 . 72833500346E+00 -0 .23405597460E-08 -0.33626434970B-10 
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